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ABSTRACT 

We present observations of the 51 Ophiuchi circumstellar disk made with 
the Keck interferometer operating in nulhng mode at N-band. We model these 
data simultaneously with VLTI-MIDI visibility data and a Spitzer IRS spectrum 
using a variety of optically-thin dust cloud models and an edge-on optically- 
thick disk model. We find that single- component optically-thin disk models and 
optically-thick disk models are inadequate to reproduce the observations, but an 
optically-thin two-component disk model can reproduce all of the major spectral 
and interferometric features. Our preferred disk model consists of an inner disk of 
blackbody grains extending to ~ 4 AU and an outer disk of small silicate grains 
extending out to ~ 1200 AU. Our model is consistent with an inner "birth" disk 
of continually colliding parent bodies producing an extended envelope of ejected 
small grains. This picture resembles the disks around Vega, AU Microscopii, and 
P Pictoris, supporting the idea that 51 Ophiuchius may be a /3 Pictoris analog. 

Subject headings: circumstellar matter — infrared: stars — interplanetary medium 
— planetary systems 
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Introduction 



^51 Ophiuchi (51 Oph), a rapidly rotating B9.5Ve star located at ISl^Jg pc ( iPerryman et al, 

19971 ). shows an infrared (IR) excess (Lm./ L^ ~ 2%) in its spectral e nergy distribution 



(SEP) due to the presence of silicate grains (IFajardo-Acosta et al.lll993l : iMeeus et al.l 12001 



Keller et al.ll2008l ). This system probably represents a rare nearby example of a young debris 
disk with gas, a planetary system just entering the late stages of formation, after the primor- 
dial gas has substantially dissipated, but before terrestrial planets have fin ished accreting. 
Its nearly edge-on disk of gas and dust and its variable absorption features (iGrady fc Silvis 
19931 : iRoberge et al.ll2002l ) suggests that 51 Oph may be an analog of P Pictoris (/? Pic). 



Several spectroscopic observations support the presence of an edge-on gaseous disk 
around 51 Oph. D ouble-peaked Ha em ission marks the presence of ionized hydrogen gas 
in Keplerian orbit fiLeinert et al.ll2004l ). Sharp CO bandhead emission features rule out a 
spherically symmetric wind and reveal a gaseous disk with tempe rature ranging from 2000 

The CO bandhead ob- 
83+^° 



- 4000 K interior to the dust sublimation radius (IThi et al. 
servations also point to a disk with inclination i 



200 



'hi et al. 



2005h or i 



( iBerthoud et al.ll2007l ). A spectral line analysis performed by lDunkin et al.l (119971 ) revealed 
a large projected rota tional velocity f o r the star o f v sin i = 267 ± 5 km s~^. Gas absorption 



features observed by iGrady fc Silvid (119931 ) and iRoberge et al.l (120021 ) are also consistent 
with an edge-on disk. 

The spati al structure of the 51 Oph dust disk remains puzzling. An HST ACS V-band 



non-d etection (jPoering et al.l 120071 ) and a Keck 18 /xm non-detection (jjayawardhana et al. 



200ll ) suggest that the disk is relatively com pact. However, far-IR photometry reveals cold 
dust grains extending out to thousands of AU ( Waters et al.lll988 : Lecavelier des Etangs et al 
1997h . 



Leinert et al.l (120041 ) obtained the first spatially- resolved measurements of the 51 Oph 



disk with the Mid-IR Interferometric (MIDI) instrument on the Very Large Telescope In- 
terferometer (VLTI). The large visibility values they measured (~ 0.65) imply that the 51 
Op h disk is relati v ely co mpact along the VLTI-MIDI projected baseline (101.2 m, 23° E of 
N). iLeinert et al.l (120041 ) fit the VLTI-MIDI visibility with a Gaussian disk and found the 
FWHM to be 7 mas, or 0.9 AU, at 12.5 /im. 

Here we present new spatially-resolved observations of 51 Oph using the Keck Interfer- 
ometer Nuller that help to c onstrain the geome try of the dust disk. We com pare these to 



the "V LTI-MIDI observations (ILeinert et al.ll2004l ) and Spitzer IRS observations (IKeller et al. 



20081 ). We simultaneously model all three data sets using two simple, edge-on disk models: 
an optically-thin model based on our zodiacal cloud and a two-layer model based on the 
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Chiang fc GoldreichI (119971 ) disk model 



Observations 



Observations of 51 Oph were conducted using the twin Keck telescopes atop Mauna 
Kea, Hawaii, operated in nulling interferometric mode on 2 Apr 2007. 51 Oph was observed 
twice in the N-band (8 - 13 /im) at an hour angle of ~ 0.5, with a projected baseline of 66.2 
m at a position angle of 47°. A calibrator star, Epsilon Ophiuchi (HD146791), was observed 
once following the target observations. Table [1] lists the details of our observations. 

The Keck Nuller operates with each K eck telescope aperture divided into two sub- 



apertures for a total of four apertures (see ISerabyn et al.l l2005l : IColavita et al.l l2008l . for 



details). Long-baseline fringes formed by combining the light from opposite telescopes are 
used to null the light from the central star and observe any spatially resolved structures. 
Short-baseline fringes formed by combining the light from two neighboring sub-apertures are 
used to remove the thermal background. 

The observable quantity is the "null leakage," the ratio of the amplitude of the short 
baseline fringe with the long baseline null fringe on target to the am plitude of the short 



baseline fringe wi th the long baseline constructive fringe on target (see lSerabyn et al.l 12005 



Barry et al.ll2008l . for details). We estimated and subtracted the systematic null leakage by 
comparing the measured null leakage of the calibration star, e Oph, with the expected null 
leakage for a limb-darkened star with the same diameter. We estimated the diameter of 
e Oph as 2.94 mas and adopted 1.5 mas error bars — much larger than the true size error 
bars — as a simple means of enforcing a systematic noise floor based on our estimate of the 
instrument performance. 



3. Data & Analysis 

3.1. Keck Nuller null leakage 

Figure [1] presents the calibrated null leakage for 51 Oph. We combined the data from 
both observations, which had nearly identical projected baselines. We limited our analyses 
to the 8 - 11.6 /im range since noise from decreased instrument throughput rendered data 
beyond 11.6 /xm unusable. For wavelengths less than 11.6 fim, the null leakage remains 
relatively fiat with a ~l-cr rise near 9.8 /im. 

We first modeled the null leakage at two different wavelengths with uniform disk and 



-4- 



Gaussian surface brightness profiles. We found angular diameters of 13.5 ± 0.5 mas and 
18.5 ± 0.4 mas at 8 and 10 fim, respectively, for the uniform disk profile. For a Gaussian 
surface brightness profile, we found FWHM of 8.1 ± 0.3 mas and 11.3 ± 0.2 mas at 8 and 
10 fim, respectively. At a distance of 131 pc, 10 mas corresponds to a transverse distance 
of 1.3 AU, suggesting that the disk is either truncated at a close circumstellar distance, or 
the axis of rotation of the near edge-on disk is oriented within a few tens of degrees of the 
projected Keck baseline. 

To better understand the geometry of the 51 Oph system, we examined our Keck Nuller 
observations together with the observations of 51 Oph made with VLTI-MIDI and Spitzer. 
Figure [2] shows a collect ion of three indepe ndent data sets from observations of 51 Oph: the 



Spitzer IRS spectrum (IKeller et al.l |2008| ) in the top panel, our N-ba nd Keck Nuller nul l 



leakage in the middle panel, and the N-band VLTI-MIDI visibility data (ILeinert et al.ll2004l ) 
in the bottom panel. 



3.2. MIDI visibility 



The VLTI-MIDI visibility was obtained on 15 Jun e 2003 with a projected baseline of 
101.2 m at a position angle of 23° ( iLeinert et al.ll2004j ). Figure [3] shows how incorporating 
this data set improves the {u, v) coverage of our analysis. Although the VLTI-MIDI baseline 
was oriented within 25° of the Keck baseline, modeled uniform disk and Gaussian surface- 
brightness profile sizes are approximately 35% smaller for VLTI-MIDI measurements than 
for Keck Nuller measurements. When we modeled the VLTI-MIDI measurements using a 
uniform surface-brightness disk, we found best -fit angular diameter s of 8.5 ±1.4 mas at 8 fim 
and 12.4 ± 1.1 mas at 10 /zm, consistent with ILeinert et al.l (120041 ). For a Gaussian model, 
we found a FWHM of 5.0 ± 0.9 mas at 8 fxm and a FWHM of 7.7 ± 0.6 mas at 10 fxm. 

The middle panel in Figure [2] shows a 1-cr rise at 9.8 fim in the Keck null leakage. The 
VLTI-MIDI visibility contains a 1-a dip at 9.7 fim. These features mirror the 10 /^m silicate 
emission feature shown in the Spitzer IRS spectrum (Section 13.31) and suggest that 51 Oph 
is more extended near 10 fim. 



3.3. Spitzer IRS spectrum 



The mid-infrared spectrum, shown in the top panel of Figure [21 was obt ained on 22 
Marc h 2004 using the Infrared Spectrograph (IRS) on the Spitzer Space Telescope ( IKeller et al, 
20081 ). Spitzer observed 51 Oph in staring mode from 5 to 36 fim using the Short-Low (SL) 
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module from 5 to 14 fim, the Short-High (SH) module from 10 to 19 /im, and the Long-High 
(LH) module from 19 to 36 fim. SL has a resolving power of R= 60-128, while SH and LH 
have R~600. 51 Oph was observed using three slit positions stepped across the nominal 
position of the sou rce. The spectral extraction and calibration methods are described in 



Keller et al.l (120081 ) . Figure [2] shows that the spectrum exhibits a pronounced 10 /im silicate 



feature and a small 18 /xm silicate feature. 

Emission from polycyclic aromatic hydrocarbon (PAH) molecules contributes to the 
mid-IR spectra of many Herbig Ae/Be st ars. We cons i dered the possibility that this emission 
could contribute to our data on 51 Oph. iKeller et al.l (120081 ) included 51 Oph in a study that 
looked for correlations between disk structure and mid-IR PAH emission from intermediate- 
mass young stars. They found no measurable emission from PAH features in the 6-13 /xm 
range in the mid-IR spectrum of 51 Oph. 



4. Modeling the 51 Oph disk 



There have been several previous attempts to model the 51 Oph IR excess. IWaters et al. 



(119881 ) fit IRAS photometric data with an optically thin, spherically symmetric dust shell 
model and found a best fit model with a d ust density proportiorial to r ~^-^ and dust tempera- 
tures ranging from ~ 100 K to ~ 1000 K. iFajardo-Acosta et al.l (|l993l ) modeled photometric 
data from the Infrared Telescope Facility (IRTF) and estimated that the IR excess could be 



attributed to astronomical 



silicates smaller than 5 /im, ranging in temperature from 400 K 



to 1000 K. iLeinert et al.l (l2004f) cornpare d an optically-thick Herbig Ae disk model with a 
puffed inner rim (iDuUemond et al.ll200ll ) to the 51 Oph MIDI visibility in the 8-14 /im 



range and IR ISO spectra and found that such a model fit poorly. 

We developed new models for 51 Oph to incorporate the new data from Spitzer and the 
Keck nuller. We do not model the detailed mineralogy of the 51 Oph disk and do not intend 
for our models to explain all of the observed spectral features to high numerical accuracy. 
We focus on the 3-D dust density distribution and disk models that qualitatively describe 
all three data sets. 

We adopted the following fitting procedure to simultaneously fit all three data sets: 



1. Model the disk contribution to the Spitzer IRS spectrum to obtain the radial distribu- 
tion of grains, the grain size, and the surface density for each model component. 



2. Using the parameters determined by the best-fit Spitzer IRS spectrum, fit the inter- 
ferometric data using 3-D models of the dust distribution to obtain the scale height of 
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each component and the disk position angle. 



To model the disk contribution to the Spitzer IRS spectrum, we first calculated the 
IR excess. We modeled the stellar contribution to the Spitzer IRS spectrum as an ideal 
bla ckbody with an effective tem perature of 10,000 K and a luminosity of 260 Lq at 131 
pc ( Ivan den Ancker et al.ll200ll ). The stellar continuum contributes on the order of a few 
percent in the N-band, so any uncertainties in the stellar luminosity or temperature are too 
small to significantly impact the interferometric or spectral responses of our models. For 
example, a 10% increase in the assumed luminosity would only raise the stellar fractional 
flux contribution from ~ 5% to ~ 7% at 10 fim. 

The 51 Oph IR excess exhibits a sharp increase in flux at wavelengths shorter than 7.5 
fim that resembles a continuum source. However, this f eature may not be con t inuum . Using 
ISO-SWS spectra with a spectral resolution of ~ .02%, Ivan den Ancker et al.l (l200ll ) showed 
this region exhibits emission features from hot circumstellar gas species, including H2O, NO, 
CO, and CO2. In light of this contamination, we chose to ignore the spectrum at wavelengths 
shorter than 7.5 fim while we were fitting our dust cloud models. We also resampled the 
Spitzer IRS data to a resolution of AA/A = 0.0185 to conserve computer time. 

To create 3-D optically-thin disk models, we used the ZODIPIC software package 



(iMoran et al.ll2004j ). This software (available online at htt p://asd.gsfc.nasa.^ov/Marc.Kuchner/home.html ) 
synthesizes images of optically thin dust disks based on the iKelsall et al.l (119981 ) empirical 
model of the solar zodiacal cloud as observed by COBE DIRBE. The grain size and inner 
and outer radius of each modeled disk component were determined by the best fit to the 
Spitzer IRS spectrum. 



4.1. Single-component optically-thin models 



For our first models we used a single component of silicate dust grains of size s distributed 
from an inner radius, Tinner, to an o uter radius, rp„tsr- We assumed a density distribution 



similar to that of the zodiacal cloud (IKelsall et al.lll998l ): the surface density S(r) oc r~°'^^. 
We numerically calculated the temperature as a function of circumstellar distance based on 
the stellar spectrum and dust optical constants, accounting for the heating of small grains 
above the local blackbody temperature. We examined 100 g rain sizes, s, r a nging from 0.05 
-2.5 fim and used the astronomical silicate emissivities from lDraine fc Led (119841 ). We used 
a non-linear least squares fitting routine to determine the best fit inner and outer disk radii 
for each grain size. 



The best-fit single-component disk model, shown in Figure HI consists of 1.0 fim grains 
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distributed from the grain sublimation radius (~ 0.65 AU) to 189 AU. As you might imagine, 
this single component does a poor job of fitting the complexities of the Spitzer IRS spectrum. 
We found that a single- component model can not adequately reproduce the relative 10 /im 
to 18 fim flux ratio, the width of the 10 fim feature, and the flux in the 13 - 15 fim range. 

4.2. Two-component optically-thin models 

An additional disk component, either of different composition, grain size, or location, 
appears necessary to fit the Spitzer IRS spectrum. Therefore, we examined several models 
consisting of two optically-thin components: an outer disk of small dust grains (< 5 /xm), 
which contribute the 10 /xm silicate emission feature, and an inner disk of large blackbody 
grains. Our optically-thin models are defined by 10 adjustable parameters: inner and outer 
radii of the outer disk (rin,i. Tout,!), inner and outer radii of the inner disk (rin,2, ''"out,2), grain 
size in the outer disk (si), scale height of the outer disk (hi/r), scale height of the inner disk 
{h2/r), surface density scaling factor of the outer disk (Si), surface density scaling factor 
of the inner disk (S2), and a common position angle for both disk components (PA). We 
assumed a fixed inner disk grain size, S2 = 100 /im, a fixed dust density radial distribution, 
n{r) oc r~^'^^ (see below), and a fixed disk inclination, i = 90°. 

First, we fit our two-component model to our resampled Spitzer IRS spectrum. We 
examined 20 values for si ranging from 0.1-3 fim and explored a wide range of inner and 
outer disk radii ranging from the dust sublimation radius for silicate grains (0.65 AU for 1 fim 
grains) to thousands of AU. Generally, for the models that best fit the spectrum, the large 
grain component stretched from ~ 0.6 AU out to ~ 5 AU and the small grain component 
stretched from a few AU out to ~ 1000 AU. 

The top panel in Figure [2] shows the total flux of the best fit model, which we refer to 
as "Model A." It also shows the contribution of each component of this model to the total 
flux. As shown in Figure [H Model A qualitatively reproduces all of the major features in 
the Spitzer IRS spectrum. 

We used the parameters that best fit the Spitzer IRS spectrum to create 3-D optically- 
thin disk models using ZODIPIC. We assumed a fixed disk inclination of 90°. We examined 
237 values of disk scale heights ranging from h/r = 0.007 to h/r = 0.2 for both the inner and 
outer disks. We examined position angles from to 180° in 1° increments. We calculated the 
Keck Nuller null leakage and VLTI-MIDI visibility for all combinations of disk scale heights 
and position angle to find the best-fit Model A, the parameters for which are listed in Table [2j 
To calculate the Keck Nuller null leakage and VLTI-MIDI visibility, we used a software suite 
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that we designed to model these instruments (see lBarry et aLll2008l ). We confirmed that the 



results from our software suite agreed with the Visibility Modeling Tool, a tool developed 



by the NExSci for simulating KIN data (http://nexsciweb.ipac.caltech.edu/vmt/vmtWeb/ ). 



The interferometric responses of Model A are shown in the middle and bottom panels 
of Figure [2l The best simultaneous fit to the Keck NuUer null leakage and VLTI-MIDI 
visibility, shown in black, illustrates that Model A does not satisfactorily reproduce both 
interferometric responses simultaneously. The Keck null leakage is well-fit, but the VLTI- 
MIDI visibility is underestimated by a factor of ~ 3, indicating that this dust disk model 
is too extended along the VLTI-MIDI baseline. Figure [2] also shows in gray the response 
corresponding to the best fit to the MIDI data alone, which does not produce a satisfactory 
Keck null leakage. 

Figure E] explains the source of the problem. The top panel in Figure E] shows the Keck 
and VLTI-MIDI 10 /im responses for Model A as a function of disk position angle. The VLTI- 
MIDI response of Model A (gray diamonds) crosses the measured VLTI-MIDI value (gray 
line) at a position angle of ~ 115°. The Keck Nuller response of Model A (black triangles) 
crosses the measured Keck null leakage (black line) at a position angle of ~ 131°. To fit 
both data sets, the 10 /im responses for the model would need to intersect their respective 
measured values at a single position angle. The top panel in this figure clearly shows that 
this does not happen; there is no single position angle for the model that works for both the 
Keck and VLTI-MIDI data. 

Figure [5] also reveals how a new model can remedy the problem. The position angles of 
the maximum in the VLTI-MIDI response and the minimum in the Keck response are fixed 
and correspond to the alignment of the instrument's projected baseline with the disk axis 
of rotation. So to improve our fit to the interferometric data, we must change the model so 
that we broaden the widths of the maximum in the VLTI-MIDI response and minimum in 
the Keck response; we must make the model more azimuthally symmetric. 

We attempted to accomplish this improvement manually by increasing the scale height 
of the disk models. However, making this change alone reduces the maximum VLTI-MIDI 
visibility of the model until it no longer reaches the measured value. So to compensate, we 
also reduced the outer radius of the inner disk; we forced the inner and outer radii of the inner 
disk to 0.54 AU and 4 AU, respectively, and re-ran the fitting procedure. We call the resulting 
best-fit "Model B." The bottom panel in Figure [5] shows that the 10 /im interferometric 
responses of Model B cross the measured values approximately simultaneously at a single 
position angle of 122°, indicated with a dotted vertical line in the figure. 

Although Model B does not fit the Spitzer IRS spectrum as well as Model A, it still 
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qualitatively reproduces the spectrum's major features, as shown in the top panel in Figure 
ini The middle and bottom panels of Figure [6] show that these changes significantly improve 
the fit to the interferometric data; Model B is consistent with both the Keck null leakage 
and VLTI-MIDI visibility. 

Figure [7] shows a simulated image of our Model B circumstellar disk at 10 fim with 1 
mas pixel resolution. The inner disk of large dust grains, truncated at 30.5 mas (4 AU), 
is the brightest feature of our model. The outer disk extends from 53 mas to ~9200 mas, 
beyond the range of the figure. 

Our models imply that 51 Oph hosts a cloud of small gr ains located at huii dreds 
to thousands of AU from t he st ar, supporting the models of IWaters et al. Jl988h and 
Lecavelier des Etangs et al. Jl997h . Our models are consistent with limits placed on the 



disk flux by previous non-detections. An HST ACS non-detection at V-band using the 1.8" 
occulting spot limits the disk surface brightness to 3 x 10~^ Jy arcsec"^ between 2" and 4" 
( iDoering et al.l 120071 ) and a Keck 18 /im nondetection places an up per limit on the surface 
brightness at 1" of 2% of the peak flux (jjayawardhana et al.ll200ll ). Our Model A, which 
extends to ~ 1200 AU, has a mean V-band surface brightness of 8 x 10^^ Jy arcsec"^ be- 
tween 2" and 4" and a mean 18 fim surface brightness at 1" of 0.2% of the peak flux, well 
within the ACS and Keck non-detection limits. Our Model B has a mean V-band surface 
brightness of 1.4 x 10~^ Jy arcsec"^ between 2" and 4" and a mean 18 fim surface brightness 
at 1" of 0.4% of the peak flux, also well within the non-detection limits. We also compared a 
10.6 /im model image against recent diffraction-limited Ke ck imaging using segment tilting 
interferometry on a single telescope (IMonnier et al.ll2009l ). Our model appears consistent 
with these data, which indicate the observed characteristic emission scale to be <30 mas, 
with >95% of emission arising within an aperture of 1.5". 

J2008h obs e rved the inner portions of the 51 Oph disk at K-band with 



TatuUi et a. 



VLTI-AMBER. iTatuUi et all fl2008f ) found the best fit to the 2.2 micron continuum visibility 
measurements using a narrow ring 0.24 AU in radius, well within the dust sublimation radius, 
and noted that this continuum may result from hot dust interior to the sublimation radius 
which is shielded frorn stellar radiation or from free-free emission from an inner gas disk 
( iMuzeroUe et al.ll2004l ). We c alculated the K-ban d visibility of our Model B for the three 
baseline orientations used by iTatuUi et al.l (120081 ) and checked these calculations with the 
Visibility Modeling Tool provided by NExScI. We calculated K-band visibilities of 0.99, 0.74, 
and 0.74 for baselines of 55 m at 34°, 82 m at 91°, and 121 m at 69°, respectively. These 
values are close to the measured visibilities of l.OiO.l, 0.8±0.05, and 0.8±0.03; the latter 
two modeled visibilities fall just below the lower limit of the measured visibilities given the 
reported uncertainties. 



- 10 - 



Because the Keck Nuller and VLTI-MIDI do not resolve these extreme inner portions 
of the disk, the K-b and observatio r is ma y probe a different structure than what the N-band 
observations probe. iTatulh et ahl (120081 ) use photometric data to estimate that the excess 
continuum contributes 25% of the total flux at 2.2 microns. By estimating the stellar flux 
at 2.2 microns and modeling the continuum flux source as a blackbody at 1500 K (dust 
at the sublimation temperature), we calculate that the continuum source's contribution at 
N-band would be no larger than 1.9 Jy. This contribution is a factor of ~5 less than the flux 
from the blackbody component of our best fit two-component optically thin disk models. 
This disparity sugg ests that an additional blackbody component is necessary beyond what 



TatuUi et al.l (120081 ) model. Additionally, our Model A and Model B blackbody components 



contribute only 1.4 Jy and 3.4 Jy o f flux at 2.2 rnicrons , respectively, less than the 3.6 Jy of 
excess continuum flux calculated by lTatuUi et al.l (120081 ). Perhaps a three-component model, 
with an outer disk of small silicate grains, an inner disk of large grains, and a continuum 
source interior to the dust sublimation radius (shielded hot dust or free-free emission from 
a gas disk), will be necessary to explain all of the interferometric and spectral data sets. 



4.3. Two-layer models 



Motivated by the lChiang &: Goldreichl (119971 ) and iDuUemond et al.l (120011 ) circumstellar 
disk models, we explored an alternative edge-on disk geometry consisting of a blackbody 
midplane layer sandwiched between two identical optically-thin surface layers of small dust 
grains. The surface layers are thin; they have a very small scale height and are unresolved 
in the direction perpendicular to the midplane. Our models were completely defined by 6 
parameters: the disk position angle, the temperature and scale height of the middle layer, 
and the temperature, grain size, and line density of the surface layers. We investigated 
100 middle-layer and 100 surface-layer temperatures ranging from 300 - 1300 K, and 6 
surface layer dust sizes from 0.1 - 1.0 /im. We used the surface-layer temperature, which 
represents the temperature of the surface layer at the outer truncation radius of the disk, to 
calculate the outer radius of the disk. Because our d i sk is e dge-on, we ignore any contribution 
by a hot inner wall. W e note that iLeinert et al.l (120041 ) fit the VLTI-MIDI data with a 
DuUemond et al.l (12001! ) flared disk model which included a hot inner wall and found this 
model to fit poorly. 

We first derived the total emitting area of each component, and therefore the line density 
of the surface layers and scale height of the middle layer, by fitting the flux of the two disk 
components to the Spitzer observations. For the 500 best fits to the Spitzer IRS spectrum, 
we calculated the Keck null leakage and VLTI-MIDI visibility as a function of position angle 
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in 1° increments. Figure [8] shows the single best fit to the three data sets. The best fit two- 
layer model qualitatively reproduces the major features of the Spitzer IRS spectrum, shown 
in the top panel of Figure El although it underestimates the fiux near 8 /im and overestimates 
the flux from 10-18 /im. The lower two panels of Figure [S] shows this model can reproduce 
the mean Keck null leakage, but cannot simultaneously reproduce the VLTI-MIDI visibility. 
The interferometric data seems to rule out this disk geometry. 



5. Discussion &z Interpretation of Models 



Our modeling efforts have yielded an optically-thin disk model that can qualitatively 
reproduce our 51 Oph observations with the Keck NuUer and also the VLTI-MIDI and Spitzer 
observations of this system. This model. Model B, is composed of an inner ring of large grains 
distributed from their sublimation radius (~ 0.5 AU) out to ~ 4 AU, and an outer ring of 
small grains (< 1 /im) distributed from an inner radius of ~7 AU to an outer radius of ~1200 
AU. We also experimented with optically-thick models, but found them unsatisfactory. 

Our preferred optically-thin disk model could be interpreted as an inner "birth" disk 
of material producing small dust grains through collisions. Grains produced in the inner 
ring that are larger than the blowout size either spiral inward under Poynting-Robertson 
drag or are collisionally fragmented. Grains smaller than the blowout size exit the system 
on hyperbolic orbits in a dynamical time; the ejected grains and marginally ejected grains 
correspond to the outer disk component. The blowout size for blackbody dust grains in a 
disk without gas is given by 

— (^ \ 

Sblowout ~ o i\/r ^ v / 

where and are the luminosity and mass of the star, res pectively, G is the g ravitational 
constant, c is the speed of light, and ps is the dust density flBurns et al.lll979l ). Assuming 
f» 260 Lq, pa 3 Mq, and Ps = 2 gm cm~^, the blowout size for 51 Oph is approximately 
50 fim, consistent with our preferred model. A similar pattern, a central population of 
large grains and an oute r ring of small grains , has be en seen in obse rvations of the debris 
disks around P Pictoris (lAugereau et al.ll200ll ). Vega (ISu et al.ll2005l ). and AU Microscopii 
fetrubbe fc Chiang! booeh . 



5.1. Gas in the 51 Oph disk 



Of course the disk around 51 Oph is not gas free (e.g. iLecavelier des Etangs et al.lll997l : 
van den Ancker et al.ll200ll : iBerthoud et al.ll2007l ). so our scenario of an inner birth ring and 
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outer disk of ejected grains requires that the gas disk is sparse enough to allow the small 
grains to be ejected from the system. Here we estimate the maximum tolerable gas density 
for which small grains are unbound in the 51 Oph system. 

The stopping time of a dust grain undergoing Epstein drag forces is given by 

tstop ~ , l^j 

PgCsound 



wher e Pg is the gas density and Csound is the sound speed (see, e.g. IWeidenschilling fc Cuzzi 



19931 ). The dynamical timescale is given by 



/ ^3 \ 1/2 



^dynamic ~ i^QM J ' ^^'^ 

where a is the semi-major axis of the grain. For the dust grains to be significantly affected 
by gas drag, tgtop < ^dynamic- Substituting Csound = i^Tg/mgY^'^, where k is the Boltzmann 
constant, Tg is the gas temperature, and rrig is the mass of an H2 molecule, we find that 
entraining a 1 /xm grain in the gas would require 

'''' ^ ^^^''T-'^'(AtiY'^'- (4) 



,1 cm-3; ~ s Vl AU. 
Assuming a gas temperature of 1400 K at a circumstellar distance of 5 AU, near the vicinity 
of the inner disk, the gas density required to stop a 1 /im grain is ~ 10^ cm~^. We conclude 
that if the disk around 51 Oph is comprised of an inner birth ring and an outer disk of 
ejected grains, then the gas density in the 51 Oph disk must be less than ~100 times the 



gas density of (3 Pic (IBrandeker et al.ll2004j ) at a circumstellar distance of ~5 AU. 



5.2. A possible warp in the disk 



The outer disk in Model B has a scale height at least ~4 times larger than the scale 
height of the inner disk. So the small grains at high orbital inclinations in the outer disk 
seem unlikely to have originated from the thin inner disk in our simple model. Such small 
grains cannot easily be perturbed to high-inclination orbits by hidden planets either, because 
of their short lifetimes compared to secular time scales. 

The difference in scale heights between the two components of our model suggests to us 
that 51 Oph may have an unobserved sub-disk of larger bodies on inclined orbits in the inner 
regions of the disk. The limited {u, v) coverage and resolution of the Keck and VLTI-MIDI 
observations could easily keep such a sub-disk hidden from our observations. This s ub-disk 
might be analogous to the X-pattern, or "warp" imaged by iGolimowski et al.l (120061 ) in the 
f3 Pictoris dis k. The large bod i es in t his sub-disk could launch the small grains onto inclined 
orbits, just as lAugereau et al.l (I2OOII ) suggested occurs in the /3 Pictoris disk. 
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5.3. Limitations of the model &; sources of confusion 

Here we discuss two possibilities that might complicate our interpretation of the mid-IR 
interferometry of 51 Oph: a) the disk composition or structure is more complex than can be 
represented by our models, or b) the disk has changed during the four year interval between 
VLTI-MIDI and Keck observations. 

The chemical composition of the 51 Oph dust disk is not well-known. Our models as- 
sumed astronomical silicates only, and ignored the possibility of a mor e complex composit ion. 



A detailed model of the 51 Oph disk composition, along the lines of iReach et al.l (120081 ) for 



example, which we leave for future work, would help to further refine our models. 

The disk geometry could also be more complex than our models can capture. As previ- 
ously discussed, there are a number of observations which suggest that an inner circumstellar 
disk is near edge-on in the 51 Oph system. Our optically-thin models of the dust disk as- 
sume an outer disk which is coplanar, but the 51 Oph disk inclination may be more complex 



Some debris disks show warps or sub-disks at different inclinations (e.g. iGolimowski et al. 



20061 ). It is possible that our assumption of coplanar disks does not allow Model A, which 
best fits the Spitzer IRS spectrum, to simultaneously fit the VLTI-MIDI visibility and Keck 
null leakage. 

The 51 Oph disk geometry may also have changed between observations. The blowout 
time for small grains originating from a belt of material at 1 AU is on the order of a single 
dynamical time (only a few years). So the outer disk could evolve significantly within the 
four years between Keck a nd VLTI-MIDI observations via a recent coUisional event (e.g. 



van den Ancker et al.l 120011 ). Additionally, 51 Oph's inner disk could feat ure complex reso- 



nant structures, such as clumps or rings due to the presence of planets (e.g. lStark fc Kuchner 



20081 ) . which rotate in and out of view on a dynamical timescale. 



A comp arison of the Spitzer IRS sp ectrum shown in Figure [2] and the ISO spectrum 



published in Ivan den Ancker et al.l (120011 ) reveals that the 5 1 Oph disk may have changed 
significantly in the 8 years between spectral observations. The silicate emission feature in 
the Spitzer IRS spectrum peaks at 10 /xm, whereas the ISO spectrum peaks at 10.5 to 11 fim. 
The slope of the right side of the silicate emission feature is also noticeably steeper in the 
Spitzer IRS spectrum; estimated slopes of the right side of the 10 /im silicate emission feature 
are ~ —1.4 Jy fim~^ and ~ —2.9 Jy /im~^ for ISO and Spitzer observations, respectively. 
These differences indicate that the 51 Oph disk may have had larger grains at the time of 
the ISO observations. Finally, the flux longward of 20 /im in the ISO spectrum is ~ 3 Jy less 
than the Spitzer IRS spectrum, possibly because there were fewer cold grains at the time of 
the ISO observations. 
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6. Summary 



We observed the 51 Oph disk at N-band using the Keck interferometer operating in 
nulhng mode. We combined the observed Keck null leakage with VLTI-MIDI visibility 
data and the Spitzer IRS spectrum and simultaneously modeled all three data sets. We 
experimented with a variety of optically-thin dust cloud models and also an edge-on optically- 
thick disk model. The Spitzer IRS spectrum ruled out the single- component optically-thin 
model, while the interferometric data ruled out our optically-thick model. 

Our preferred model consists of two separate populations of large and small grains. The 
three data sets are best simultaneously fit by our Model B (Table E]). This model, shown in 
Figure [6|, contains a disk of larger grains that extends from the grain sublimation radius out 
to ~ 4 AU and a disk of 0.1 fim grains that extends from ~ 7 AU to ~ 1200 AU. 

This model may be consistent with an inner "birth" disk of continually colliding parent 
bodies. The large grains (> 50 fim) produced by the parent bodies make up the inner disk, 
while the small grains (< 50 fim) are blown outward and eventually ejected from the system 
by radiation pressure. The large scale height of the outer disk compared to the inner disk 
suggests that the small grains which compose the outer disk may originate from an unseen 
inclined sub-disk or from a population of inclined cometary bodies. 

Although the 51 Oph disk seems puzzling at first, perhaps it is not so strange after 
all. The distribution of grain sizes i n our models is not un ique to the 51 Oph sy stem, but 
has been ob served in the Pictoris (lAugereau et al.ll200lf ). Vega fISu et al.ll2005f ). and AU 
Microscopii (IStrubbe fc Chiang |2006| ) disks. Our models suggest that the 51 Oph disk may 
be another member of a class of debris disks which exhibit similar dust distributions. Our 
models also indicate that there may be two sources of dust at different inchnations around 
51 Oph. These models, together with prev ious observations of variable absorption features 
( 1 Grady fc Silvid Il993l : iRoberge et al.l |2002| ) which may be due to transient infalling bodies 
suggest that the 51 Oph dust disk may well be an example of a /? Pictoris-like system. 
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Table 1. Keck NuUer Observation Log 



Object 


Type 


Time 


U 


V 


Air Mass 






(UT) 


(m) 


(m) 




51 Oph 


target 


15;08;24 


52.00 


48.97 


1.39 


51 Oph 


target 


15;09;26 


51.90 


48.87 


1.39 


e Oph 


caUbrator 


15:37:21 


37.34 


59.76 


1.21 



Table 2. Best-fit optically-thin model parameters with 99.73% joint confidence estimates 





(/im) 


nn,i 
(AU) 


''out,! 

(AU) 


hi/r 


Sl 
(Zodis)'^ 


S2 


nn,2 
(AU) 


»'out,2 

(AU) 


h2/r 


(Zodis)"^ 


PA 

n 


Model A 


n 1 +0.06 


2 44+0-07 


-,200+2^00 


0.03t 


1.15iro|xlo5 


100^ 


U. (1D_Q QQg 


14 5+°-'' 


U.UOU_Q 


2.49lg 0^x105 


101+0.15 
-^■^-^-0.05 


Model B 


O.lt 


7 1+0.3 
'■^-0.3 


lZUU_gQQ 


0.19* 


2.o5i;^:[j^xio5 


100^ 


0.54'' 


4.0^ 


0.04t 


4.3t°:Jxl05 


1 99+0.5 
J-^^_0.15 



Tixcd parameter in both models 
''Fixed value in Model B 

'^One "zodi" refers to a face-on optical depth of 10"^ at 1 AU 
t Upper limit 
* Lower limit 
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Fig. 1. — Keck NuUer null leakage measurements for 51 Oph. 



- 20 - 




Wavelength (fim) 



Best siimiltaneous fit 
Best fit to MIDI data 



/h 2^^^ tsK- 



Wavelenglh (|im) 



Best simultaneous fit 
Best fit to MIDI data 




Wavelength (|im) 



Fig. 2. — Model A, the best fit to the three data sets using a 2-component optically thin 
disk model. Top panel: the observed Spitzer IRS spectrum (solid black line), and best-fit 
modeled spectrum. Middle panel: the observed Keck NuUer null leakage (triangles) shown 
with error bars, the best-fit modeled null leakage when simultaneously fitting both Keck and 
VLTI-MIDI data sets (black line), and the best-fit null leakage when fitting the VLTI-MIDI 
data set alone (gray line). Bottom panel: the observed VLTI-MIDI visibility (diamonds) 
shown with error bars and modeled VLTI-MIDI visibilities. The parameters of Model A 
which best fit the two interferometric data sets simultaneously are listed in Table [2l The 
best fit to the VLTI-MIDI data alone is shown in gray and corresponds to a disk with a 
position angle of 117°, an outer disk scale height of hi/r ^ 0.038, and an inner disk scale 
height of h2/r ^ 0.024. 
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Fig. 3. — {u, v) coverage of the observations discussed in this paper. The outer boundary of 
the Spitzer {u, v) coverage is shown as a circle, corresponding to the Spitzer mirror diameter 
at 8 microns. 
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Fig. 4. — Best fit single-component disk model to the Spitzer spectrum (sofid black line). 
The best fit features 1 jivn grains, distributed from 0.65 AU to 189 AU. 
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Fig. 5. — Top panel: Modeled interferometric response at 10 //m as a function of position 
angle for the Keck NuUer (black triangles) and VLTI-MIDI (gray diamonds) for Model A. 
The solid black and gray lines show the measured 10 /xm responses for Keck and VLTI- 
MIDI, respectively, along with shaded regions representing their respective measurement 
uncertainties. The modeled Keck and VLTI-MIDI responses do not cross their measured 
values together at any one position angle. Bottom panel: Modeled interferometric response 
at 10 fjLm for Model B. The modeled responses cross the corresponding measured values at 
approximately 122°. 



-24- 




Fig. 6. — Model B, the best fit to the three data sets using a 2-component optically thin 
disk model where the inner and outer radii of the inner disk are forced to 0.54 AU and 4 
AU, respectively. The best-fit parameters of Model B are listed in Table El 
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Fig. 7. — Simulated 10 /xm image of Model B in units of flux per pixel, with a pixel size of 1 
mas. The disk midplane is oriented at 122° East of North. The parameters of Model B are 
listed in Table [2l 
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Fig. 8. — Best simultaneous fit to the three data sets using a two-layer optically-thick disk 
model. The best fit model is a disk oriented with a position angle of 38° and truncated at 
outer radius 3.4 AU. The disk features a middle layer of 360 K blackbody grains with a disk 
height of 0.163 AU from the midplane, and 2 surface layers of 0.7 iim grains heated to 835 
K, each with line densities of 8.2 x 10^^ AU"^ 



